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ABSTRACT 

We observe and characterize the scattering of acoustic wave packets by a sunspot, in a regime 
where the wavelength is comparable to the size of the sunspot. Spatial maps of wave traveltimes and 
amplitudes are measured from the cross-covariance function of the random wave field. The averaging 
procedure is such that incoming wave packets are plane wave packets. Observations show that the 
magnitude of the traveltime perturbation caused by the sunspot diminishes as waves propagate away 
from the sunspot - a finite-wavelength phenomenon known as wavefront healing. Observations also 
show a reduction of the amplitude of the waves after their passage through the sunspot. A significant 
fraction of this amplitude reduction is due to the defocusing of wave energy by the fast wave-speed 
perturbation introduced by the sunspot. This "geometrical attenuation" will contribute to the wave 
amplitude reduction in addition to the physical absorption of waves. In addition, we observe an 
enhancement of wave amplitude away from the central path: diffracted rays intersect with unperturbed 
rays (caustics) and wavefronts fold and triplicate. Thus we find that ray tracing is useful to interpret 
these phenomena, although it cannot explain wavefront healing. 

Subject headings: Sun: helioseismology — Sun: interior — Sun: oscillations — Sun: photosphere — 
sunspots — Sun: surface magnetism 



1. INTRODUCTION 

The propagation of solar seismic waves is affected by 
sunspots. The scattering phase shifts and wave absorp- 
tion coefficients can be measured by several techniques 
such as Fourier-Hankel analysis and time-d istance analy- 
sis (see recent review bv lGizon et al.l l2010 and references 
t herein). 

ICameron. Gizon. fc Duval (l200l (hereinafter CGD) 
studied the interaction of /-mode plane-wave packets 
with a sunspot using time-distance helioseismology. As 
we shall describe briefly below, the data averaging strat- 
egy used by CGD greatly reduces the noise and enables 
a detailed study of the waveforms far from the scatter- 
ing region. In this study, we extend the observations 
to modes with different radial orders and wavelengths 
that are closer to the sunspot radi us, in order to investi- 
gate finite- wavele ngth effects (e.g. iNolet fe Dahlenl [20001 : 
iHung et all 120011 ). Not only do we measure traveltime 
shifts but also wave amplitude perturbations with re- 
spect to the quiet Sun. 

We use one-minute cadence Doppler velocity im- 
ages measured by the Michelson Doppler Imager 
(jScherrer e t al. 1995) onboard the Solar and Heliospheric 
Observatory (SOHO/MDI). The images are remapped 
using Postel's azimuthal equidistant projection in order 
to trac k the motion of the sunspot in Active R egion 9787 
(CGD: lGlion et al.ll2lM iMoradi et alllMlCft . The map 
scale is 0.12 heliographic degrees per pixel. Nine data 
sets are obtained, one for each day of the period 20-28 
January 2002. The selected sunspot is nearly circular and 
does not evolve over the duration of the observations; its 
outer-penumbral radius is R — 20 Mm. 

Filters are applied to the Dopplergrams in 3D Fourier 
space to extract waves with the same radial order (ridge 
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filters). Five different filters are applied, for modes / 
through p4 (corresponding to radial orders n — 0-4). 
Table [T] lists the characteristics of the wave packets for 
each radial order, including the central frequency, v, and 
the central wavelength, A. Broadly speaking, the finite- 
wavelength effects are expected to be observable in the 
regime where R/X approaches unity (this condition will 
be refined in Section [4|. These effects should thus be 
easier to see with the wave packet. 

2. CROSS-COVARIANCE FUNCTION AND WAVE PACKETS 

The cross-covariance of the filtered Doppler velocity is 
computed in the same way as in CGD: 



C n (x,y,t)= / (f> n (t')ct>n(x,y,t' + t)dt', 



(1) 



where (x, y) is a local Cartesian coordinate system with 
origin at the center of the sunspot, t is the time, <p n is 
the filtered Doppler velocity for modes with radial or- 
der n, and 4> n is the average of 0„ over the line L at 
x = —43.73 Mm. The spacial averaging over L is equiv- 
alent to filtering out all the wavevectors that are not 
perpendicular to L: only waves that propagate toward 
and away from L contribute to </>„. 
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Figure 1. Cross-correlation function F^(x, y,T) (solid line) be- 
tween the cross-covariance C4,(x,y,t) at (x,y) = (60,0) Mm and 
the corresponding quiet-Sun cross-covariance C? . The envelope 
(dashed line) is determined by demodulation. The traveltime shift, 
St, and the wave packet amplitude, A, are indicated by arrows. In 
practice, we fit the peak nearest to r = with a parabola to mea- 
sure the time shift. In this particular example, the time shift and 
reduced amplitude are St = —57 s and A = 0.44 with respect to 
quiet Sun. 

As explained bv iGizon et all (|2010l ). for t > 0, the 
cross-covariance is used to follow the propagation of two 
plane wave packets initially at L and propagating to- 
ward +x and —x respectively. The interaction of the 
waves with the sunspot will occur in the (x > 0,t > 0) 
quadrant. 

To reduce noise, the cross-covariances are averaged 
over nine days and over azimuthal angle about the 
sunspot center (see CGD). A n example pl o t of c ross- 
covariances has been given by iMoradi et all poll Fig- 
ure 17). 

3. WAVE TRAVELTIMES AND AMPLITUDES 

We define the traveltime shift, Sr(x,y), as the time- 
lag that maximizes the similarity between the mea- 
sured cross-covariance, C n (x,y,t), and a sliding quiet- 
Sun cross-covariance, C% s (x,t). More precisely, 5t is the 
time r that maximizes the function 

P , ^ _ JC n (x,y,t)C^(x,t-r)dt 

*»V>v> T >- f\c?{x,t)\*at ■ [) 

The quiet-Sun cross-covariance C% s (x,t) is constructed 
by averaging C n (x,y,t) over \y\ > 100 Mm, far away 
from the region of wave scattering, and it is zeroed out for 
t < 0. The cross-correlation tr aveltime shift defin ed here 
is also used in geophysics (e. g. iZaroli et al.ll2010D and is 
analogous to the definition of lGizon fc Birchl ((2002). The 
cross-correlation F n (x, y, r) at each spatial point (x, y) as 
a function of r is demodulated using the Hilbert trans- 
form to obtain its envelope. The maximum, A{x, y), of 
the envelope of F n (x, y , r) is a measure of the wave packet 
amplitude. By definition, A is equal to one in the quiet 
Sun. An example measurement of St and A for C4 at 
(x, y) — (60, 0) Mm (after interaction with the sunspot) 
is shown in Figure [T] 

For each wavepacket, we obtain spatial maps of the 
traveltime shift Sr(x,y) and wave amplitude A{x,y), 
shown in Figure [2] The maps of 5t clearly show nega- 
tive traveltime anomalies behind the sunspot along y = 



while the maps of A show the amplitude reductions. 

We only display the results for distances from L larger 
than three wavelengths in each panel in order to exclude 
the near field where the t < and t > branches of C n 
are not well separated. This is not a limitation of this 
work as we are interested in wave scattering in the far 
field. For / and pi modes there are positive traveltime 
anomalies on the left side of the penumbra, which are 
likely to be caused by the outward-directed flows in the 
penumbra and the moat. 

4. WAVEFRONT HEALING AND FINITE- WAVELENGTH 
EFFECTS 

The left panel of Figure [3] shows the traveltime shifts 
as a function of x along y = 0. For all the modes, 
the magnitude of the traveltime shift gradually de- 
creases as waves propagate away from t he sunspot. This 
phenomenon, called wa vefront healing (|Nolet fc Dahlenl 
l2000t iHung et al.l l2001). is due to the decreasing ampli- 
tude of the scattered wavefield with distance. This is a 
finite-wavelength effect that cannot be explained by ray 
theory, which w ould give a consta nt traveltime perturba- 
tion for x > R (IGizon et al.ll2006[) . 

The occurrence of finite-wavelength effects at a mea- 
surement point depends on whether the diameter of 
the sunspot is smaller than the width of the first Fres- 
nel zone of that point. For an observation point at 
(x, y = 0) (the point that enters the computation of 
the cross-correlation), the width (measured along the 
y-axis) of the first plane-wave Fresnel zone is given by 
Ayp(x) — 2y/(x + A /4) 2 — x*, assuming a path length 
perturbation of A/4 (|Gudmundssonl lT996). Let us define 
xf, the characteristic x such that the sunspot fills the 
first Fresnel-zone width, i.e. Ayp(xp) — 2R. For each 
wave packet, the value of xf is provided in Table [T] The 
corresponding first plane-wave Fresnel zone is bounded 
by a parabola with focus at point (xf, 0) (white curves in 
top panels of Figure [5]). In the regime x > xp, the first 
Fresnel zone is always larger than the sunspot: finite- 
wavelength effects dominate and the ray approximation 
is no longer valid. 

The right panel of Figure [3] shows p4-mode traveltime 
shifts as a function of y, at four different values of x. 
As the distance away from the sunspot increases, we 
see wavefront healing around y = 0. In addition, we 
see that the width of the travel-time perturbation in the 
transverse direction grows with x; this observation can- 
not be explained in the context of linearized ray theory, 
according to which travel-time perturbations are com- 
puted along unpert urbed paths (y = const.) by using 
Fermat's principle (H ung et al.ll200il ). This is another 
warning against using linearized ray theory in sunspot 
seismology. 

As we shall see in the next section, a proper ray tracing 
calculation indicates that a fast wave-speed anomaly (like 
a sunspot) would defocus the rays and lead to travel-time 
perturbations that spread in the transverse direction as x 
increases. However, ray tracing will still overestimate the 
traveltime shifts in the far field and would be incapable 
of explaining wavefront healing along y — 0. 

5. AMPLITUDE ENHANCEMENTS AND CAUSTICS 

The bottom panels of Figure [2] show that the transmit- 
ted wave packets have a reduced amplitude around y = 
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Figure 2. Spatial maps of traveltime shifts 8t(x, y) (top) and wave packet amplitudes A(x, y) (bottom) with respect to quiet Sun. From 
left to right are results for wavepackets /, pi, p2, P3, and P4. Waves are initially located at L (red lines) and propagate toward +x. The 
black circles denote the outer boundary of the penumbra. The parabolic curves (white lines) depict the boundaries of the first plane-wave 
Fresnel zones of a point located at (x, y) = (xp, 0), chosen such that the sunspot fills the width of the Fresnel zone. We do not display the 
results in the near fields (gray regions where distances from L are smaller than 3A). 
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Figure 3. Observed traveltime shifts, St. The left panel shows the traveltime shifts along y = for modes / through p4, which heal 
as the distance from the sunspot, x, increases. The right panel shows the traveltime shifts of p4 modes as functions of y at four different 
values of x. Notice the spread in the transverse direction as x increases. The hatched and cross-hatched area indicate the locations of the 
penumbra and umbra respectively. Data are averaged over patches of size 16 X 16 Mm 2 and error bars give the standard deviation of the 



compared to the quiet-Sun value, for all radial orders. 
It is known that wave absorption by sun spots causes 
a reduction in outgoing wave amplitude (jBraun et al.l 
119871 119881 ). as the result of partial mode conversion of 
incoming waves into slow magn etoacoustic waves that 
propagate down the su nspot (e.g. lSpruit &: Bogdanlll992t 
ICallv fc Bogdan1[T997l) . 

In our case (plane wave geometry), a significant frac- 
tion of the amplitude reduction that is observed on-axis 
behind the sunspot is due to the defocusing of wave en- 
ergy by the fast wave-speed perturbation introduced by 
the sunspot. This can be illustrated by a simple 2D 
ray tracing experiment shown in the left panel of Fig- 
ure HJ This "geometrical attenuation" will contribute to 
the wave amplitude reduction in addition to the physi- 
cal absorption of waves. Their respective contributions 
would require better modeling as ray tracing cannot take 
account of wavefront healing. 

Apart from the on-axis amplitude reduction, an en- 



hancement of wave amplitude is seen away from the axis 
y = for modes pi 1 P3, and p^ (bottom panels in Fig- 
ure [5]) . This off-path enhancement is more pronounced 
for the p4 wave packet (right panel of Figure 3} . Like 
above, we suggest that this off-path feature is caused by 
sunspot-induced refraction. Since rays tend to bend to- 
ward the slow-speed medium, diffracted rays will cross 
over the unperturbed rays on both sides off the central 
axis. The envelope of the bundle of intersected rays is 
called a caustic in optics. The wavefronts fold and trip- 
licate as they pass through the caustics. Inside the caus- 
tics, all three arrivals will contribute to cross-correlation 
traveltime measurements and the interpretation of trav- 
eltim e becomes ambiguous under the ray-the oretical pic- 
ture (jNolet fc DahlenlhoOOtlHung et al.»200TI ). The loca- 
tion of the caustics depends on the perturbation strength 
and the sunspot geometry. In the case of p4 modes, a 10% 
wave-speed perturbation provides a rather satisfying ex- 
planation for the on-axis power deficit and off-path en- 
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Figure 4. The wave amplitude perturbations for p4 caused by the sunspot. Left: Map of observed amplitudes A (color scale). The 
white lines arc 2D rays traced through a toy sunspot model, such that the wave speed is enhanced by 10% inside the circle and transitions 
smoothly to the background value (u\ = 53.28 km s" 1 ) at the dashed circle. The wavefronts are given by the black lines. The thick black 
lines indicate the boundaries of the causics where rays intersect and wavefronts fold. Right: Observed amplitudes A as a function of y at 
fixed distances from the sunspot ranging from x = 50 to 200 Mm. Notice the amplitude enhancements around y ~ ±60 Mm. 



hancement (left panel of Figured]). We note that the off- 
path enhancements are also seen in the cross-covariance 
power maps (not shown here), though not as clearly. 



6. DISCUSSION 

We have measured the traveltime shifts and amplitudes 
of plane wave packets as they traverse a sunspot using 
a cross-covariance technique. Observations show on-axis 
wavefront healing and amplitude reduction. Traveltime 
anomalies spread in the transverse direction with dis- 
tance and amplitude enhancements are seen away from 
the central ray. 

Wavefront healing is a finite-wavelength phenomenon. 
The only other direct observational evidence of finite- 
wavelength effects in tim e -dista nce helioseismology was 
reported by iDuvall et al.l (|2006ft who studied the inter- 
action of waves with sub-wavelength magnetic features. 
Wavefront healing cannot be modelled by ray theory. 
This implies that ray-based interpretations are not al- 
ways appropriate in sunspot seismology. 

We have seen that ray-tracing (unlike linearized ray 
theory) is quite useful to study the geometrical atten- 
uation due to the defocusing of wave energy by a fast 
wave-speed anomaly. We conclude that wave amplitude 
reduction is not necessarily due to wave absorption by 
the sunspot. Ray tracing is also useful to interpret off- 
path amplitude enhancements due to wavefront folding, 
although it cannot explain wavefront healing. 

The observations of traveltime shifts and wave am- 
plitude that have been presented in this paper contain 
a wealth of information about the seismic signature of 
sunspots. Their interpretation requires methods of anal- 
ysis more sophisticated than linearized ray theory (which 
is still often used in sunspot seismology today). 
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